Orbital polarons in the metal-insulator transition of manganites 
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The metal-insulator transition in manganites is strongly influenced by the concentration of holes 
present in the system. Based upon an orbitally degenerate Mott-Hubbard model we analyze two 
possible localization scenarios to account for this doping dependence: First, we rule out that the 
transition is initiated by a disorder-order crossover in the orbital sector, showing that its effect on 
charge itineracy is only small. Second, we introduce the idea of orbital polarons originating from 
a strong polarization of orbitals in the vicinity of holes. Considering this direct coupling between 
charge and orbital degree of freedom in addition to lattice effects we are able to explain well the 
phase diagram of manganites for low and intermediate hole concentrations. 

PACS number(s): 72.80.Ga, 71.30.+h, 71.38.+i, 75.50.Cc 



I. INTRODUCTION 

The doping dependence of the properties of manganese 
oxides poses some of the most interesting open problems 
in the physics of these compounds. First to be noticed 
is the peculiar asymmetry of the phase diagram that is 
most pronounced in the charge sector: Regions of high 
(x > 0.5) and low (x < 0.5) concentration of holes are 
characterized by such contrasting phenomena as charge 
ordering and metalicity, respectively. In the latter re- 
gion — which we wish to focus on — the metallic state 
can be turned into an insulating one by raising the tem- 
perature above the Curie temperature Tc- Introducing 
the notion of double exchange which associates the rel- 
ative orientation of localized Mn ^g spins with the mo- 
bility of itinerant e g electrons, early work has identified 
this transition to be controlled by the loss of ferromag- 
netic order inherent to the metallic state.lil It is believed 
that lattice effects are also of crucial importance in this 
transition. Within the lattice-polaron double-exchange 
picturep'0 the crossover from metallic to insulating be- 
havior is controlled by the ratio of polaron binding energy 
Eb to the kinetic energy £"kin of charge carriers: 



A 



E b 

•Ekin 



(1) 



When forming a bound state with the lattice, charge car- 
riers loose part of their kinetic energy. Hence, polarons 
are stable only if this loss in energy is more than com- 
pensated by the gain in binding energy, i.e., if A > 1. In 
a double-exchange system, this critical coupling strength 
may be reached by raising the temperature — the double- 
exchange mechanism then acts to reduce the kinetic en- 
ergy and hence to increase A. Spin disorder and spin- 



polar on e, 
above T^ J 



ects further enhance the carrier localization 
The doping dependence of the metal-insulator 



transition, however, is not fully captured in this picture. 



Namely the complete breakdown of metalicity at hole 
concentrations below a; cr i t ~ 0.15 - 0.2 that occurs despite 
the fact that ferromagnctism is fully sustained remains 
an open problem which we address in this paper. 

The effective coupling constant A in Eq. (fil) has origi- 
nally been introduced for non-interacting electrons. The 
itinerant e g electrons in manganites, on the other hand, 
are subject to strong on-site repulsions which necessitates 
to accommodate the definition of A. According to numer- 
ical studiesjj the basic physics underlying Eq. (fil) remains 
valid even in correlated systems: As in the free-electron 
case, the metal-insulator transition is controlled by the 
competition between the polaron binding energy and the 
kinetic energy of charge carriers. Nevertheless, correla- 
tion effects might renormalize these two relevant energy 
scales, presumably introducing a doping dependence. In 
fact, the Gutzwiller bandwidth of correlated electrons 
scales with the concentration of doped holes; one could 
therefore be inclined to set E\^ n oc xt, where t denotes the 
hopping amplitude. But this approach reaches too short: 
The Gutzwiller picture describes only the average kinetic 
energy of the system. In contrast, the relevant quan- 
tity for localizing the holes doped into a Mott-Hubbard 
system is the characteristic saergy scale of charge fluc- 
tuations which remains ex tu Pictorially this quantity 
corresponds to the kinetic energy of a single hole. We 
thus conclude that a more thorough treatment of corre- 
lation effects is needed in order to explain the peculiar 
doping dependence of the metal-insulator transition in 
manganites. 

In this paper we analyze two mechanisms that could 
drive the localization of charge carriers at small hole con- 
centrations x. First we explore the possibility of the 
metal-insulator transition to be controlled by a disorder- 
order crossover in the e g -orbital sector. The idea is the 
following: Orbital fluctuations are induced by the mo- 
tion of holes and hence possess an energy scale oc xt. 



At large x, orbitals fluctuate strongly and inter-site or- 
bital correlations are weak. As the concentration of holes 
is reduced, fluctuations slow down until a critical value 
of x has been reached — promoted by Jahn- Teller and 
superexchange coupling, an orbital-lattice ordered state 
now evolves. We analyze the extent to which this tran- 
sition in the orbital-lattice sector affects the itineracy of 
holes. Finding almost similar values for the kinetic en- 
ergy of doped holes in orbitally ordered and disordered 
states we are lead to conclude that the development of 
orbital-lattice order is not sufficient to trigger the local- 
ization process. Next we turn to analyze a second sce- 
nario of the metal-insulator transition for which we in- 
troduce the concept of orbital polarons. Similar to spin 
polarons in correlated spin systems, orbital polarons are 
a natural consequence of strong electron correlations and 
the double degeneracy of on-site levels — in mangan- 
ites the latter follows from the degeneracy of e g orbitals. 
We argue that holes polarize the orbital state of e g elec- 
trons on neighboring sites: A splitting of orbital levels is 
evoked by a displacement of oxygen ions and also by the 
Coulomb force exerted by the positively charged holes. 
Being comparable in magnitude to the kinetic energy of 
holes, the orbital-hole binding energy can be large enough 
for holes and surrounding orbitals to form a bound state. 
The important point is that the stability of these orbital 
polarons competes not only against the kinetic energy of 
holes but also against the fluctuation rate ex xt of or- 
bitals: The faster the latter fluctuate, the less favorable 
it is to form a bound state in which orbitals have to give 
up part of their fluctuation energy. Combining this new 
orbital-polaron picture with that of conventional lattice 
polarons we are able to explain well the phase diagram 
of manganites at low and intermediate doping levels. 



II. ORBITAL DISORDER-ORDER TRANSITION 

In this section we analyze the impact of a disorder- 
order transition in the orbital-lattice sector onto the itin- 
eracy of holes. Our motivation is that a sudden freezing 
out of orbital fluctuations below a critical doping concen- 
tration could significantly impede the motion of holes, 
hence initiating the metal-insulator transition. By com- 
paring the bandwidth of holes both in orbitally ordered 
and disordered states, we are able to refute this idea: The 
orbital sector is shown to have only little influence onto 
the charge mobility in manganites. 



A. Disordered State 



The model accounts for the presence of strong on-site re- 
pulsions that forbids more than one e g electron to occupy 
the same Mn site as well as for the double-degeneracy of 
e g levels. At low temperatures and intermediate doping 
levels, the double-exchange mechanism induces a paral- 
lel alignment of spins. Treating deviations from this fer- 
romagnetic ground state only on a mean-field level as 
is discussed below, the core spins can be discarded and 
the spin indices of e g electrons may be dropped; the t-J 
Hamiltonian then becomes&tj 
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with z = 6. Nearest-neighbor bonds along spatial direc- 
tions 7 G {x,y, z\ are denoted by (ij) 1 - We use con- 
strained operators c ia — c ia (1 — ni) which create an e g 
electron at site i in orbital a only under the condition 
that the site is empty. The first term in Eq. (g) de- 
scribes the inter-site transfer of constrained e g electrons. 
The transfer amplitude depends upon the orientation of 
orbitals at a given bond as is reflected by the transfer 
matrices 
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a representation with respect to the orbital basis a G 



{|3; 



y 2 )} has been chosen here. Due to 



its non-diagonal structure, orbital quantum numbers are 
not conserved by Hamiltonian (0) - inter-site transfer 
processes induce fluctuations in the orbital sector. The 
second term in Eq. (0) accounts for processes involv- 
ing the virtual occupation of sites by two e g electrons. 
This superexchange mechanism establishes an inter-site 
coupling between orbital pseudospins of overall strength 
J = zt 2 /U, where U is the on-site repulsion between 
spin-parallel e g electrons. The pseudospin operators are 



r x /y 



](<r z ±V3a x ), 



1 



(3) 



with Pauli matrices a x J z acting on the orbital subspace. 
Jahn- Teller phonons mediate an additional interaction 
between orbital pseudospins which is of the exact same 
form as the superexchange term. The numerical value of 
J has to be chosen such as to comprise both effects. We fi- 
nally note that deviations from the ferromagnetic ground 
state underlying Hamiltonian (0) are treated within con- 
ventional double-exchange theory. The transfer ampli- 
tudfj-i then depends on the normalized magnetization m 



We begin by investigating the bandwidth of holes in 
a strongly fluctuating, orbitally disordered state. Our 
starting point is the t-J model of double-degenerate e g 
electrons which, via Hund's coupling, interact ferromag- 
netically with an array of localized S = 3/2 core spins. 



t = t Qy /{l + m 2 )/2, 



(4) 



where to denotes the hopping amplitude between spin- 
parallel Mn sites. 

To observe the strongly correlated nature of e g elec- 
trons, it is convenient to introduce separate particles 



for charge and orbital degrees of freedom.113 The metal- 
lic phase of manganites can be well described within an 
orbital-liquid picture that accounts for orbital fluctua- 
tions by emploMiag a slave-boson representation of elec- 
tron operators SO 

C ia = fia^i- 

Here orbital pseudospin is carried by fermionic orbitons 
fi a and charge by bosonic holons bi. Introducing mean- 
field parameters \ = t^ 1 Y, a /3 t^ULfip) and x = ( & ! & i)» 
where x is the concentration of holes and x ~ \ > the two 
types of quasiparticles can be decoupled: 

/U, -U+^)E'f (/1/^ + H.c), (5) 
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(ij) 



(6) 



Diagonalizing the above expressions in the momentum 
representation one obtains 
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with index v = ± and dispersion functions 



^ = (xt 



uj k = 6%i 
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1 



^eo(fe) 



where eo(fe) = c a 
e 2 (fc) = V3(c x - 



£l (fc) 



,)/2 



c y )/2 with c 7 = cosfc 7 . The essence 



of this slaved-particle mean-field treatment is that or- 
bital and charge fluctuations are assigned different energy 
scales. This is reflected by the bandwidths of orbiton and 
holon quasiparticles, respectively: 



W orb 



6xt - 
6i. 



J, 



(7) 



The former quantity W or b sets the energy scale of orbital 
fluctuations — the terms proportional to xt and J de- 
scribe fluctuations induced by the motion of holes and 
by the coupling between pseudospins, respectively. The 
latter quantity Whi finally defines the itineracy of holes 
in the orbital-liquid state. The variation of the holon 
bandwidth with the onset of orbital order is in the focus 
of our interest in the remainder of this section. 



B. Instability Toward Orbital Order 



coupling between orbitals on neighboring sites which in- 
troduces a bias towards orbital-lattice ordering. Compet- 
ing against the energy scale of orbital fluctuations oc xt, 
order in the orbital sector is expected to evolve below a 
critical doping concentration ir cr it. We investigate this 
instability of the orbital-liquid state by introducing the 
inter-site coupling term 



z *-^ 



e e iQ~, 



(9) 



with z = 6 and Tj = (sin Oaf + cos 0cr?)/2 acting on the 
orbital subspace. We note that Eq. (|)|) is a simplifica- 
tion of the superexchange coupling term in Hamiltonian 
(|2J) — internal frustration makes the latter difficult to 
handle. For = tt/2 and Q = (tt,tt,0), the pseudospin 
interaction in Eq. (0) favors a staggered-type orbital ori- 
entation 



\e) ± = (\3z 2 -^)±\x 2 -y 2 ))/V2 



(10) 



within x-y planes repeating itself along the z direction; 
this closely resembleSpthe type of order observed experi- 
mentally in LaMnOsJlj The breakdown of the orbitally 
disordered state, i.e., the development of orbital order, 
is signaled by a singularity in the static orbital suscepti- 
bility (cQcrf.g) u;= o. Employing a random-phase approx- 
imation, the latter can be expressed as 



(°* 
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(11) 



with vertex function J q = J (cos q x + cos q y — cos q z ) /3 
and the shorthand notation ((J x (t x )q = (^q<J : 1q)uj=o- 
Bare orbital susceptibilities (. . .)° are evaluated using or- 
biton propagators associated with the mean-field Hamil- 
tonian (||). Numerically solving for the pole in Eq. (|ll|), 
we find the following expression for the critical doping 
concentration: 



J_ 



(12) 



which is valid for a; cr it <C 0.5. At concentrations below 
this critical value an orbitally ordered state is to be ex- 
pected. With J = 0.13 eV as estimated from the struc- 
tural phase transition observed in LaMnC>3 at T = 780 K 
(Ref. [l2[ ) and t — 0.36 eV (Ref. ||) we obtain a critical 
doping concentration of x cr ;t = 9%. This result indicates 
that the metallic state of manganites is indeed instable 
towards orbital-lattice ordering at doping concentrations 
that are not too far from those at which the system is 
observed to become insulating. 



The above treatment of orbital and charge fluctuations 
is based upon the notion of a strongly fluctuating orbital 
state that is far from any instability towards orbital or- 
der. In real systems such instabilities do exist: Jahn- 
Teller phonons and superexchange processes mediate a 



C. Ordered State 

Up to this point we have studied the bandwidth of 
holes in an orbitally disordered state and the instability 




State 1 



State 2 




FIG. 1. Coherent (a) and incoherent (b) hole motion in 
antiferro-type orbital order. Incoherent processes involve the 
creation of an orbital excitation of energy J. 



of the system towards orbital-lattice order. In the follow- 
ing we analyze to which extent the itineracy of holes is af- 
fected by this disorder-order transition in the orbital sec- 
tor. Namely we are interested in the bandwidth of holes 
moving through an orbitally ordered state; this quantity 
is then compared to our previous result of Eq. m) for an 
orbital liquid. Foremost an important difference between 
models with orbital pseudospin and conventional spin is 
to be noticed here: In the latter systems, spin is con- 
served when electrons hop between sites. This implies 
that hole motion is constrained in a staggered spin back- 
ground. In contrast, the transfer Hamiltonian (pi) of the 
orbital model is non-diagonal in orbital pseudospin — an 
orbital basis in which all three transfer matrices i"' 3 , t" 13 , 
and if%r are of diagonal structure does not exist. This 
allows holes to move coherently even within an antiferro- 
type orbital arrangement [see Fig. 11(a)]. For this reason 
only a moderate suppression of the hole bandwidth is to 
be expected in the presence of orbital order. We calcu- 
late the bandwidth for the specific type of orbital order 
introduced in Eq. ( J10[ ) . Starting from the transfer part of 
Hamiltonian (0) and keeping only the hopping matrix el- 
ements that allow for a coherent movement of holes (i.e., 
projecting out all orbitals which do not comply with the 
ordered state) we obtain 



IT, 



coh 
hi 



At. 



(13) 



This result indicates a reduction of the holon bandwidth 
by w 30% as compared to the disordered state [see Eq. 
(0)]. While not being dramatic, a quenching of the band- 
width by one third should nevertheless be sufficient to 
induce the localization process, e.g., via the formation of 
lattice polarons. However, it is important to note that 
Eq. (13) accounts solely for the coherent motion of holes; 
incoherent processes involving the creation and subse- 
quent absorption of orbital excitations are neglected [see 
Fig. [j](b)]. Only if the ordered state is robust, i.e., if it 
costs a large amount of energy for an electron to occupy 
an orbital that does not comply with the long-range or- 



bital alignment, these incoherent processes become negli- 
gible. This limit does not apply to manganites where the 
orbital excitation energy is J < t only. Before a conclu- 
sion about the role of orbital order in the metal-insulator 
transition can be drawn, these incoherent processes have 
to be investigated in more detail. 

In the following we study the influence of incoherent 
processes onto the motion of holes, employing an "orbital 
wave" approximation: Starting from the assumption that 
long-range orbital order has developed and that fluctua- 
tions around this ordered state are weak, we use a slave- 
fcrmion representation of the electron operators in the 
transfer Hamiltonian (0): 

A — ftt f 

Within this picture, the orbital pseudospin is assigned 
to bosonic orbitons and charge to fermionic holons. The 
lattice is then divided into two sublattices which are as- 
cribed different preferred pseudospin directions [see Eq. 
©]: 



T =(|3z 2 - 
4^(|3z 2 - 



\ x ~ V )) /v2 on sublattice A, 



If 



'■)) /V2 on sublattice B. 



In analogy to conventional spin-wave theory, excitations 
around this ground state can be treated by employing the 
following mapping of orbiton operators bi a onto "orbital- 
wave" operators [3f. 



3iT 



1 sublattice A, 
Pi sublattice B, 



5 a 



sublattice A, 
sublattice B. 



In the momentum representation the transfer Hamilto- 
nian (0) then becomes 

H t = Y J Uhftfh + E M + 7fc+p/3- P ] ftfh+p. (14) 



kp 

-t(c x + c y - 2c,j)/2 and -y k = t[(2 - \/3)c x + 
f 2c z ]/2 with c 7 = cos fc 7 . The first term in 
describes the coherent motion of holes within a 
At. The second term describes 
the interaction of holes with excitations of the orbital 
background. The dynamics of the latter is controlled 
by the inter-site coupling Hamiltonian (|9|) which in the 
momentum representation becomes 



Here ujk = 

(2 + V3)c y 

Eq. (|£ 

band of width IT™' 1 



H, 



J 



\PtPk 



(15) 



Hamiltonian ( |15[) describes dispersionless, non- 
propagating orbital excitations of energy J. The local 
nature of orbital excitations follows from the absence of 
frustration effects in the inter-site orbital coupling term 

The interaction of holes with orbital degrees of freedom 
changes the character of the hole motion: Scattering on 
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FIG. 2. Born approximation of the holon self energy: The 
double line denotes the "dressed" holon propagator with 
self-energy contributions, the dashed line represents the "or- 
bital-wave" Green's function. 



orbital excitations leads to a suppression of the coherent 
quasiparticle weight and a simultaneous widening_n£_the 
holon band. In analogy to studies of spin systems,BEJ we 
analyze these effects by employing a self-consistent Born 
approximation for the self energy of holes — within this 
method all non-crossing diagrams of the self energy are 
summed up to infinite order while crossing diagrams are 
neglected. Restricting ourselves to the case of a single 
hole moving at the bottom of the band, we obtain the 
following expression for the holon self energy (see Fig. 

1): 



Yi{iui) 



:t 2 ^- 2 



^7pG(lUJ-J, P ). 
P 



(16) 



The Matsubara frequencies are defined as iu = i(2n + 
1)ttT, where T denotes temperature and n an integer 
number. 

Our first aim is to study the loss of coherency in the 
hole motion. This caa.be done by employing a dominant- 
pole approximation:^ We split the holon propagator in 
Eq. (Hq) into its coherent and incoherent parts, 



G{iu,k) 
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(17) 



where a^ denotes the quasiparticle weight and u>k the 
not-yet-known renormalized holon dispersion. Keep- 
ing only the coherent part and using a*. = [1 — 
{d / du)Y,' {u)}~ 1 with S'(a;) = Re[£(w -> u + i0+)], the 
following recursion relation for the quasiparticle weight 
is obtained: 



-i -i 



i + t 2 Tif 



p P (^-w p - Jf 
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with q = (0,0, 7t) at the bottom of the band. Equation 
( |l8| ) can be approximately solved by expanding the inte- 
grand around q, which yields 




for J > t, 



(19) 



for J < £. 



In the limit J/t — ► oo, orbitals become static and co- 
herent hole motion with a q — 1 is recovered. In the 
opposite limit J/t — > 0, the holon quasiparticle weight is 



0.2 



3, 
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FIG. 3. Spectral function of a hole moving at the bottom 
of the band. Different values of J/t are used: The spectrum 
is completely incoherent for J — 0, having a lower bound at 
w m in = — 3t. With increasing values of J/t, spectral weight is 
shifted from the incoherent part of the spectrum to a coherent 
quasiparticle peak (denoted by a vertical line). In the limit 
J/t — ► oo, the quasiparticle peak is at ojqp = —2t and has 
accumulated all spectral weight. 



completely lost, indicating strong scattering of holes on 
orbital fluctuations. 

Next we turn to study the renormalization of the holon 
bandwidth. Inserting G(iu>, k) = [iu — Wfc — 'E(iu))]^ 1 into 
Eq. (Eq) leads to a recursion relation for the holon self 
energy: 



Y,(iu) = t 2 



"iv 



p 



J — aYi(iu — J) 



(20) 



The factor a = (z — I)/ z partially accounts for the con- 
straint that forbids more than one orbital excitation per 
site — the hole may therefore not return to a previously 
visited site unless to reabsorb an excitation. We solve 
Eq. d2Q) numerically and determine the spectral function 
Pq(u) — — ~Im[G(iu — > u + i0 + , q)] at the bottom of the 
band. The result is shown in Fig. g. Different values of 
J/t are used. In the limit J/t — > 0, the spectrum is com- 
pletely incoherent and extends down to w m i n = — 3t cott£~ 
sponding to a holon bandwidth of Wu = 2|w m i n | = 6t&3 
In this limit the hole creates its own disorder and ef- 
fectively moves within an orbital-liquid state character- 
ized by strong orbital fluctuations. In the opposite case 



J/t — > oo, all spectral weight accumulates in a quasi- 
particle peak (denoted by a vertical line) at wqp = — It 
which corresponds to a bandwidth of Whi = 2|wqp| = At. 
The orbital state is static here and excitations are com- 
pletely suppressed. At finite values of J/t, the total spec- 
tral weight is divided into coherent and incoherent parts. 
The latter is separated from the quasiparticle peak by the 
orbital excitation energy J. Processes in which the hole 
creates more than one orbital excitation are reflected by 
a succession of peaks in the incoherent spectrum. For 
J = t/3 which is realistic to manganites, the quasiparti- 
cle peak accounts for w 65% of the spectral weight and 
the width of the holon band is 




Whi « 5.7*. 



(21) 



FIG. 4. Polarization of e g levels on sites next to a hole: 
Phonons and Coulomb interaction induce a splitting of energy 

A = A ph + A ch . 



Comparing the above number with our previous result 
Whi = 6i for the orbital-liquid state, we find a reduction 
of about 5%. We therefore conclude that a disorder- 
order crossover in the orbital sector has only a secondary 
effect on the kinetic energy of charge carriers, ruling it 
out as a possible driving mechanism to initiate the metal- 
insulator transition in manganites. 



III. ORBITAL POLARONS 

In the preceeding section we have considered charge 
carriers to interact with the orbital sector via the transfer 
part of Hamiltonian (J2j). While this coupling was shown 
to be responsible for a shift of spectral weight from the 
coherent to the incoherent part of the holon spectrum, 
the effect onto the full bandwidth was found to be only 
small. In the following we point out that in an orbitally 
degenerate Mott-Hubbard system there also exists a di- 
rect coupling between holes and orbitals stemming from 
a polarization of e g orbitals in the neighborhood of a 
hole. This coupling is strong enough for holes to form a 
bound state with the surrounding orbitals at low doping 
concentrations. Based upon this picture we show that a 
strong reduction of the bandwidth comes into effect as 
orbital-hole bound states begin to form. 



A. Polarization of Orbitals 

The cubic symmetry of perovskite manganites is lo- 
cally broken in the vicinity of holes which results in a 
lifting of the e g degeneracy on sites adjacent to a hole 
(see Fig. [|). Here we discuss two mechanisms that are 
foremost responsible for this level splitting: (1) a dis- 
placement of oxygen ions that move towards the empty 
site; and (2) the Stark splitting of e g states which is in- 
duced by the Coulomb force between "positive" hole and 
negative electrons. The magnitude of the degeneracy lift- 
ing A = A ph + A ch is estimated as follows: The former 
phonon contribution A ph originates in the coupling of 



holes to the lattice breathing mode Q\ and of e g elec- 
trons to two Jahn- Teller modes Q2 and Q3: 

#ei-ph = - Yl {diQun'l + gzQtiOi + gaQaiof 



K 



-Qll 



(22) 



where n' 1 denotes the number operator for holes and 

the Pauli matrices erf act on the orbital subspace. 
The coupling constants are g\ and gi rs 53 and K 
is the lattice spring constant. Hamiltonian (G2) medi- 
ates an interaction between empty and occupied sites. 
The effective Hamiltonian describing this coupling is ob- 
tained by integrating Eq. (g2|) over oxygen displacements 
Qi = (Qui Q21, Qsi)- For a given bond along the z direc- 
tion this yields H z = -|A ph v^o* with 



A ph = gi g 2 V2/(3K) « (gi/ga)EjT. 



(23) 



A lower bound for this quantity is giAjea by the Jahn- 
Teller energy, i.e., A ph > Ej T « 0.2 eV]B assuming that 
coupling to the breathing mode is at least as strong as 
coupling to the Jahn- Teller modes. Next we estimate the 
contribution to the e g -level splitting that follows from the 
Coulomb interaction between a positively charged hole 
and an e g electron on a neighboring site. The magnitude 
A ch of this splitting is assessed by taking into account 
the covalence of Mn 3d and O 2p orbitals, which gives 



A 



ch 



3 2 R 

t7 %ii-Mii- 



(24) 



The covalency factor 7 = t pc i/A. p d can be obtained from 
the transfer amplitude and the charge gap between Mn 
and O sites, t pd w 1.8 eV and A pd « 4.5 eV.fLlI Together 
with a lattice spacing of i?Mn-Mn = 3.9 A this leads to 
A w 0.4 eV. In total the polarization of e g levels on sites 
next to a hole yields an energy splitting A rs 0.6 eV. Be- 
ing comparable in magnitude to the transfer amplitude t 
this number strongly indicates a direct coupling of charge 
and orbital degrees of freedom to be of importance in 
manganites. 




FIG. 5. Orbital polaron in the strong-coupling limit: Six 
e g states point towards a central hole. 



The splitting of e g levels effects all six sites surrounding 
a hole. From the above Hamiltonian for a bond along 
the z axis, analogous expressions for x and y directions 
are derived by a rotation in orbital isospin space. The 
complete Hamiltonian for the cubic system is then 



ffch-orb = -A ^ n> i T j, 



(25) 



with orbital pseudospin operators given by Eq. ||). 
Hamiltonian (Eq) promotes the formation of orbital po- 
larons. For low enough hole concentrations these con- 
sist of a bound state between a central hole and the 
surrounding e g orbitals pointing towards the hole as is 
shown in Fig. pi. This configuration also yields a large 
amplitude of virtual excursions of e g electrons onto the 
empty site. Thus, besides minimizing the interaction en- 
ergy of Hamiltonian (£5|), it also allows to lower the ki- 
netic energy. We note that these virtual hopping pro- 
cesses locally enhance the magnetic moments of core and 
e g spins via the double-exchange mechanism, providing 
a large effective spin of the orbital polaron. This natu- 
rally explains the development of ferromagnetic clusters 
experimentally observed at temperatures above Tc-ej 



B. Binding Energy 

In conventional lattice-polaron theory, the binding en- 
ergy is a function of the coupling constant g and the 
stiffness of the lattice which is controlled by the spring 
constant K: The energy gain stemming from the in- 
teraction between charge carriers and the lattice com- 
petes against the deformation energy of the crystal. In 
the case of orbital polarons, the underlying picture is 
very similar. Here the coupling constant is given by 
the orbital-charge interaction energy A, while the energy 
scale Worb = Qxt + J is a measure of the "stiffness" of 
the orbital sector. These two quantities are expected to 
determine the binding energy of the orbital polaron: 



I 



• 



(b) 



xt 



I i I 



FIG. 6. (a) Orbital fluctuations with an energy scale oc xt 
and (b) inter-site correlations oc J compete against the or- 
bital-hole binding energy A. 

The role of W£ rh can be illustrated as follows (see Fig. 
|6j): In an orbital- liquid state, orbitals have to give up 
part of their fluctuation energy in order to form a bound 
state with a hole. As a consequence, polarons are stable 
only if orbital fluctuations are weak. Furthermore, po- 
larons have a frustrating effect on inter-site orbital cor- 
relations. The local orientation of orbitals favored by 
Hamiltonian fl2q) does in general not comply with the 
orientation that would minimize the Jahn- Teller and su- 
perexchange coupling between orbitals on nearest- and 
next-nearest-neighbor sites of the hole. Thus, in order 
to form a bound state orbitals have to give up part of 
their inter-site correlation energy J as well. The fact 
that the polaron binding energy is controlled by the or- 
bital energy scale W OI b — Qxt + J has direct implications 
for the phase diagram of manganites: Due to the dop- 
ing dependence of W rb, orbital polarons are stable only 
at low hole concentrations where fluctuations are weak. 
The tendency of the system to form polarons is therefore 
most pronounced in the lower part of the phase diagram. 
To derive an expression for the polaron binding energy, 
the following approach is used: First we consider a static 
hole placed in an ordered state without fluctuations. We 
then calculate the reduction of the total energy due to the 
interaction Hamiltonian (Eq). All approximations made 
in the following aim at discarding terms reminiscent of a 
specific type of orbital order, while preserving the most 
general structure of the orbital-hole binding energy. We 
focus on a single site located next to the hole in the z 
direction. The orbital configuration at this site is de- 
termined by the coupling to the hole described by Eq. 
( pqj as well as by the orientation of neighboring orbitals 
which couple via superexchange and Jahn- Teller effect; 
the latter interaction is determined by Hamiltonian (M). 
Treating all orbitals except the one explicitly considered 
here on a mean-field level, the following Hamiltonian is 
obtained for the selected site: 



HI = - (At| 



Jrf) 



(27) 



E° b rh = /(A, Worb). 



(26) 



Here rf = (sin Oaf + cos0of)/2 fixes the orbital orien- 
tation which would minimize the interaction energy with 
the orbital background — for the type of order used in 
Eq. (J10|), e.g., 9 = ±7r/2. In general, this orientation 
does not coincide with the \3z 2 — r 2 ) configuration fa- 
vored by the orbital polaron which is described by the 
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FIG. 7. Polaron binding energy E% T : The solid curve is a 
plot of Eq. ( |2q ) for the orbitally ordered state as a function of 
J (upper axis), the dashed curve corresponds to Eq. (B4) for 
the orbital-liquid state and is plotted as a function of Worb 
(lower axis). 



first term in brackets. The state actually chosen by the 
system then depends on the ratio of A and J as well as on 
the angle O. We determine the energy of this state and, 
since we are not interested in any specific type of order, 
average over the angle 8. Finally subtracting the A — » 
limit and multiplying with the number of bonds connect- 
ing the hole to its surrounding, the following expression 
for the polaron binding energy is obtained: 



£ 6 orb = 3 [\/A 2 + ,P - J 



(28) 



In the limit J/A — > oo, the orbital state is very stiff and 
cannot be polarized by the hole; the binding energy then 
vanishes as 3A 2 /(2J). On the other hand, the polariza- 
tion is complete in the limit J/A — > 0, yielding a max- 
imum value E£ rh = 3A for the binding energy. A plot 
of E™ h is shown in Fig. [?] (upper curve). We note that 
the functional form of Eq. ( pf ) differs from the conven- 
tional lattice-polaron case where E% = g 2 /(2K). This 
is due to the fact that there exists an upper limit of the 
orbital polarization in which the orbitals around a hole 
have been fully reoriented to point towards the empty site 
(see Fig. H); technically the existence of this upper bound 
is reflected by the hard-core nature of the Pauli opera- 
tors in Eq. (|27f . The familiar form of the binding energy 
oc A 2 J J is recovered only for the weak-coupling limit in 
which orbital distortions around the hole are small. 

Next we consider a static hole placed in a strongly fluc- 
tuating orbital-liquid state. Hamiltonian (Eq) imposes a 



splitting of e g levels on the sites next to the hole. The 
orbiton quasiparticles of Sec. 01 scatter on these local po- 
tentials, which may lead to the formation of an orbiton- 



hole bound state. To calculate the polaron stabilization 
energy we again consider a single site next to the hole 
in the z direction. The local potential imposed by the 
close-by hole is of the form 



m = -a-. 



6 ■ 



(29) 



where if = \o\. We calculate the effect of successive 
scattering of orbitons on the above potential employing a 
T-matrix formalism; interference between different scat- 
tering centers is neglected here. The correction to the 
orbiton Green's function that seizes the effect of Hamil- 
tonian fl29| ) is given by 

8G olh (w; R, R) = G Q olb (iLu; R, 8)T s (iuj)G OTh (iu; 8, R), 



(30) 



with the scattering matrix 



Ts(iu>) 



cr 2 A/2 



l-cr*AG° rb M/2- 



Here G® Tb (iu>; R, R) denotes the orbiton propagator 
of the system in the absence of the scattering poten- 
tial. The elements of this 2x2 matrix are given by 
lG° orh (iu;;RR)} a P = -(TV/W^L and are con- 
trolled by the mean-field Hamiltonian (]5|). The on-site 
Green's function is G° rb M =J^r[G° lb (ico;S,8)]. Inte- 
grating over lattice sites, Eq. (J30|) becomes 



JO(iu) = -f 



1 



l-AG°(iu>)/2 l + AG°(ito)/2 



, (31) 



with SG OT b(iu>) — 2 J2r Tr[(5G or b(iw; R, R)]. The change 
in the total energy of the system which is induced by the 
scattering potential can now be obtained from Eq. ( |3l| ) 
by employing 



SE 



dw (to — /i) Sp(uj), 



(32) 



where Sp(oj) = —(l/ir)Im[5G(iuj — > u> + i0 + )] denotes 
the scattering contribution to the density of states. The 
orbiton chemical potential is set to /i = in the follow- 
ing. We evaluate Eq. (32) by approximating the on-site 
Green's function by 



G°M 



1 



W 



■In 



iu) - W orb /2 



iu + Worb/2 



(33) 



This expression yields a constant density of states for 
the translationally invariant system which resembles the 
result that can be obtained numerically from the mean- 
field Hamiltonian (ra) . Approximately solving the integral 
in Eq. (32) and multiplying the result with the number 
of nearest neighbors of the hole, we finally arrive at the 
following expression for the polaron binding energy: 



^orb 



3A 



y(cihy- 1) + 



yln2 



y 



(tt/2)* 



(34) 



with y = Worb/A. Equation (|34j) describes the polaron 
stabilization energy in a strongly fluctuating orbital- 
liquid state. A plot of this function is shown in Fig. 
ffl (lower curve). The binding energy reaches its maxi- 
mum E£ lh — 3A if orbital fluctuations are weak, while it 
vanishes as 31n2(A 2 /W rb) in the opposite limit where 
fluctuations are strong. 

In the last two paragraphs we have calculated the 
polaron binding energy in an orbitally ordered as well 
as in an orbital-liquid state. Although very different 
approaches were used to describe these complementary 
cases, the expressions obtained closely coincide: As can 
be seen in Fig. ffl, the two functions are almost indistin- 
guishable if one identifies J in Eq. (|2q ) with vW rb in Eq. 
(p4), where v ss 0.72 is a numerical fitting factor. Moti- 
vated by the observation that W or b = 6xt + J reduces to 
J in the limit of strong orbital correlations J 3> xt, we 
discard the fitting factor in the following by setting v = 1 . 
We believe these two cases to be smoothly connected as 
should come out in a more elaborate treatment of the 
problem. Based upon these considerations, we conclude 
that Eq. ( p8| ) can be used to model the polaron bind- 
ing energy of both the fluctuating and the static orbital 
state: 



K h 



3(1 



A 2 



wL - w 



(35) 



where W ov h = 6xt + J in former and W OI b = J in the lat- 
ter case. We finally note that in deriving expressions ( J2q ) 
and (pj) for the binding energy, all six sites surrounding 
the static hole were assumed to be occupied. Since at 
finite doping the probability of a site being occupied is 
only (1 — x), we renormalize Eq. ( pq ) by this average 
occupation factor. 

To summarize, E™ h in Eq. ( p5| ) represents the energy 
to be gained by polarizing the orbital background around 
a static hole. This number depends on the orbital energy 
scale W or b = Qxt + J, i.e., on orbital fluctuations and 
inter-site orbital correlations which both tend to suppress 
the polaron binding energy. W or h is to be considered as 
the counterpart of the lattice stiffness K in conventional 
polaron theory. The important difference is that W or b ex- 
plicitly depends on x which has important consequences 
for the phase diagram of manganites: Orbital polarons 
can form only at low doping concentrations where orbital 
fluctuations are weak and the binding energy is conse- 
quently large. 



C. Polaron Bandwidth 

The orbitally degenerate Mott-Hubbard system is in- 
stable towards the formation of orbital-hole bound states 
at low doping concentrations. In this small-polaron 



regime, holes are pinned by the binding potential and can 
move only if being thermally activated. At low temper- 
atures these processes can be neglected; coherent charge 
motion is then possible solely due to quantum tunneling. 
Since the polaron is a composite object consisting of a 
hole and several orbitals, the amplitude of these tunnel- 
ing processes is expected to be small as is indeed shown 
in this section. 

Here we restrict ourselves to polarons moving in an 
orbitally ordered state. Our analysis is based upon the 
following idea: By allowing a hole to polarize the sur- 
rounding orbitals, the system reduces its ground state 
energy by E£ lh . This energy is mostly lost if the hole 
hops to another site since the orbital sector cannot im- 
mediately adopt to the new position — orbital fluctua- 
tions are slow compared to those of holes. After a short 
time the system returns to the ground state, most likely 
by transferring the hole back to its original location. But 
there is also a small probability for the hole to keep its 
new position while the orbital sector adapts to the re- 
location. This is possible due to the non-orthogonality 
of configurations in which orbitals point towards the old 
and the new location of the hole, respectively. The po- 
laron tunneling amplitude is then given by the transfer 
amplitude t of holes multiplied by the overlap between 
states with orbitals pointing towards the old and new 
position of the hole, respectively. This overlap is calcu- 
lated as follows: We use Hamiltonian ( p7|) to determine 
the orientation of a single orbital next to the hole - 
all other orbitals are treated on a mean-field level. As 
the hole hops, this orbital has to change its orientation 
from pointing towards the hole to being aligned with the 
background. The projection between these two states is 



P 



1 

V2 



1 + 



,/ 



VA 2 + J 2 



1/2 



(36) 



An average over the angle O specifying the type of orbital 
order in Eq. (E7J) has been performed here. Other orbitals 
undergo the reversed process: Originally being aligned 
with the background, they turn towards the hole as the 
latter hops onto a neighboring site. In total there are 
2(z -l)»;2z orbitals that have to reorient. The overlap 
between the initial and the final state is then given by 
P = p 2z , yielding 



P 



We rewrite Eq. 



P = 



1 



./ 



VA 2 + J 2 



(37) 



1 - 



exp 



VA 2 + J 2 - J 




2VA 2 + J 2 




VA 2 + J 2 - 


J 


VA 2 + J 2 





(38) 



where the exponential form becomes exact for large co- 
ordination number which has been set to z = 6. The 
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FIG. 8. Ratio of polaron and holon bandwidth e~ v . The 
solid line is based on the exponential form in Eq. (p9|), the 
dashed line represents Eq. (fcSTJ). The polaron bandwidth 
shrinks which decreasing values of J/ A associated with an 
enhanced orbital polarizability. 



denominator in the exponent is identified as the orbital 
binding energy given by Eq. (Eq). Hence we finally arrive 
at the following compact expression for P = e~ v : 



-n — 



exp 



E orb 



VA 2 + J 2 



(39) 



A plot of this function is shown in Fig. |8|. The physi- 
cal significance of e~ v is that it relates the holon to the 
polaron bandwidth: 



Wp 0l = Win e" 



(40) 



As the system becomes critical towards the formation 
of orbital-hole bound states, polarons replace holes as 
charge carriers. Our result shows that this transition is 
accompanied by an exponential suppression of the band- 
width. Strictly speaking the translationally invariant sys- 
tem remains a metal; in reality, however, the small band- 
width makes polarons susceptible to localization, e.g., by 
trapping in the random potential of impurities. The sup- 
pression of the bandwidth is most pronounced if the po- 
laron binding energy is large: The orbitals around the 
hole are then strongly distorted, which necessitates a sig- 
nificant reorientation to allow the hole to hop. We note 
that the expression in Eq. (p9|) is similar to the result 
obtained in conventional latticc-polaron theory where 



3 q = exp 



-E, 



ph , 



v 6 /u>o). Here Eg* = g 2 /(2K) denotes 

the polaron binding energy and ujq the phonon frequency; 
the latter corresponds to (A 2 + J 2 ) 1 / 2 in our orbital- 
polaron theory. 

Equation (|39|) has been derived for a static, non- 
fluctuating orbital state. Following the discussion in de- 
riving the polaron binding energy, we generalize the re- 
sult to account for orbital fluctuations as well. This is 
done by replacing the inter-site correlation energy J by 



the more general orbital energy scale W or h 
Hence we obtain 



exp 



£ orb 



v/A 2 + W 2 



6xt + J. 



(41) 



where the polaron binding energy E£ rb is now given by 
Eq. (H). 

To summarize, the development of orbital polarons 
leads to a sharp reduction of the bandwidth. In this 
regime the orbital-hole bound state can move only as an 
entity via quantum-tunneling processes. Since the po- 
laron extends over several lattice sites, the transfer am- 
plitude is exponentially small. The bandwidth reduc- 
tion is controlled by the ratio E^' h /[A 2 + W 2 rh } 1/2 which 
is a measure of the orbital distortions around a hole. 
Strong orbital fluctuations and inter-site orbital corre- 
lations weaken the polaron effect by suppressing these 
distortions. 



IV. METAL-INSULATOR TRANSITION 

As was shown in the preceeding section, the forma- 
tion of orbital-hole bound states leads to an exponential 
suppression of the bandwidth which makes the system 
prone to localization. In a double-exchange system, this 
crossover from a free-carrier to a small-polaron picture 
can be initiated either by a reduction of the doping con- 
centration or by an increase in temperature; the former 
acts via an enhancement of the polaron binding energy, 
the latter by constricting the motion of holes via the 
double-exchange mechanism. In this section we combine 
our orbital-polaron picture with the theory of conven- 
tional lattice polarons to develop a scheme of the metal- 
insulator transition in manganites. 

The transition from a free-carrier to a small-polaron 
picture is governed by the dimensionlcss coupling con- 
stant A rb = E° lh /Dhi where E^ h is the polaron binding 



energy given by Eq. (|35|) and Z?hi = VFhi/2 = 3t is the 
half-bandwidth of holes. The lattice breathing mode of 
Eq. ( p2|) adds an additional contribution A p h = E% /-Dhi 
with E^ = g 2 /(2K), which further promotes the forma- 
tion of polarons. The coupling constant thus becomes 



A 



E orb 



E, 



pli 



D 



hi 



(42) 



The critical value that separates free-carrier and small- 
polaron regimes is A = 1. For A » 1 small polarons 
have fully developed and the bandwidth is reduced by an 
exponential factor 



exp 



jgorb 



VA 2 + W 2 




(43) 



with 7 = 1. We note that interference effects between 
orbital and lattice coupling are neglected in Eq. (f43). 
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FIG. 9. Crossover from the free-carrier to 

the small-polarons regime. At A < 1 polarons are unstable 
and the bare holon bandwidth is recovered. At A > 1 orbital 
polarons form, leading to an exponential suppression of the 
bandwidth. To avoid an abrupt drop of the bandwidth at 
A = 1, Eq. (J4J) is used to connect the two regimes. 

For A < 1 the free-carrier picture is recovered and holes 
move in a band of width Whi = 6i. This implies 7 = 
in Eq. rt43), yielding e~ v = 1. To simulate the crossover 
between the two regimes, we phenomcnologically employ 
the function 

InWl + ffl, 

A 2 (44) 

for A < 1, 

with (3 = [1 — 1/A 2 ] 1 / 2 . This function has been proposed 
for strongly coupled electron-phonon systems (see, e.g., 
Ref. n9) and avoids an unphysical sudden drop of the 
bandwidth as A = 1 is reached. The crossover hence 
obtained for e~ v is depicted in Fig. Q. 

Up to this point we have mainly focused on the role of 
the orbital energy scale W or b in the formation of small 
polarons. We now turn to analyze in more detail the ef- 
fect of temperature. The latter controls the bandwidth 
Whi via the double-exchange mechanism: At low tem- 
peratures all spins are ferromagnetically aligned and the 
transfer amplitude reaches its maximum. With increas- 
ing temperature the ferromagnetic moment weakens, con- 
stricting the motion of charge carriers. Specifically, the 
transfer amplitude changes with the normalized magne- 
tization m as 



i = W(1 + ™ 2 )A 



(45) 



The magnetization depends on temperature via the self- 
consistent equation 



Bg(am) with a 



35 T c 



1 + 5 T ' 



(46) 
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denotes the Brillouin function. The average magnetic 
moment per site varies with doping: The moments 5 C = 
3/2 and s — 1/2 of core and e g spins combine on average 
as 



s = | + i(i-*). 



(47) 



Finally, the Curie temperature Tc in Eq. (|46|) is con- 
trolled by the strength J c s of ferromagnetic exchange 
bonds via 



Tc^^S(S- 



l)Jeff. 



(48) 



The fitting parameter v compensates for an overestima- 
tion of Tq in the mean-field treatment. Double-exchange 
as well as superexchange processes are responsible for 
establishing the ferromagnetic links between sites. The 
magnitude|-Q|f this coupling in the limit of large Hund's 
coupling istll 



Jefi = 



1 
I25 2 



x(l - x)xte~ n + (1 - x) 



■ 2x 2 4 
u 



(49) 



where 



The first term in squared brackets of Eq. ( ff9[ ) stems 
from the coherent motion of holes/polarons and repre- 
sents the conventional double-exchange contribution to 
Tc- The factor e~ v accounts for the rescaling of the 
coherent bandwidth as the small-polaron regime is en- 
tered. The second term is due to superexchange pro- 
cesses. It describes the high-energy virtual hopping of e g 
electrons which is insensitive to a polaronic reduction of 
the bandwidth. It is also noticed that superexchange in 
an orbitally degenerate system is of ferromagnetic na- 
ture because-iif _the large Hund's coupling present in 
manganitesOlHj'EJ Superexchange hence dominates the 
ferromagnetic interaction in the small-polaron regime. 

The system of equations presented above controls the 
electronic and magnetic behavior of manganites at low 
and intermediate doping levels. A critical coupling A = 1 
leading to the formation of polarons can be reached ei- 
ther by lowering the doping concentration or by increas- 
ing the temperature — the former enhances E% , the 
latter quenches Whi- The equations are interrelated and 
have to be solved recursively. As a result of this self 
consistency, the breakdown of the metallic bandwidth at 
A = 1 is expected to be rather sharp: With the evolu- 
tion of small polarons, the coherent band width shrinks, 
thereby weakening the magnetic exchange links. Double 
exchange then drives the system even farther towards the 
strong-coupling limit. 



V. COMPARISON WITH EXPERIMENT 

To illustrate the interplay between the system of equa- 
tions presented in the proceeding section, we numerically 
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FIG. 10. Magnetic and electronic phase diagram obtained 
within the present theory. Dashed and dot-dashed lines rep- 
resent magnetic and electronic transitions, respectively, the 
simultaneous transition in both channels is denoted by a solid 
line. The phases are: paramagnetic insulator (PI), paramag- 
netic metal (PM), ferromagnetic metal (FM), and ferromag- 
netic insulator (Fl). 



extract from them the T-x phase diagram. While it is 
obvious that T = Tc is a suitable criterion to sepa- 
rate the low-temperature ferromagnetic from the high- 
temperature paramagnetic state, more care has to be 
taken to distinguish between metallic and insulating be- 
havior. Our theory describes the reduction of the band- 
width which follows from the formation of small polarons. 
However, strictly speaking, the system remains metallic 
even in the strong-coupling limit since polarons can still 
move by tunneling. It is therefore necessary to define a 
critical value of the bandwidth beyond which additional 
effects such as pinning to impurities are implicitly as- 
sumed to set in and finally turn the system into an insu- 
lator. The specific criterion used here is only of marginal 
importance, as feedback effects discussed above induce a 
quick collapse of the bandwidth once a critical coupling 
A = 1 is reached. For simplicity we define A < 1 to be a 
metal and A > 1 to be an insulator. 

The following parameters are chosen for calculating the 
phase diagram: The orbital polarization energy is set to 
A = 0.55 eV, yielding a binding energy comparable to 
the phononic one E% — 0.45 eV; the phonon frequency 
is ojq = 0.05 eV, the interaction between orbitals_J = 
0.13 eV, the bare transfer amplitude to = 0.36 eV,H and 
U = 4.0 eV. The fitting parameter v = 0.55 (Ref. ||) 
is adjusted to reproduce the values of Tc observed for 
Lai-^Sr^MnOs.u The result is shown in Fig. [10| 

Our most important observation is that the doping de- 
pendence of orbital polarons makes the system more in- 
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FIG. 11. Variation of the normalized magnetization m with 
temperature. At low hole concentrations (upper curve), the 
system is close to a transition to the small-polaron regime — 
a small increase in temperature then already sufficiently im- 
pedes the hole motion for polarons to form, resulting in a 
sharp drop of the magnetization. At large x (lower curve), 
the conventional double-exchange picture is recovered. 



sulating at low and more metallic at high doping levels. 
Convincingly this is seen in the complete absence of met- 
alicity at x < 0.15 and the appearance of a metallic phase 
above Tc at x > 0.4. The region 0.15 < x < 0.4 in which 
colossal magnetoresistance is experimentally observed is 
characterized by a simultaneous magnetic and electronic 
transition-irom a ferromagnetic metal to a paramagnetic 
insulator O The role of polarons in this transition is 
most pronounced at low hole concentrations. This can 
be seen from the behavior of the magnetization as Tc 
is approached from below (see Fig. O): At low x, the 
polaron binding energy is large and the system is close 
to localization. A small reduction of the bandwidth via 
double exchange is then sufficient to trigger the forma- 
tion of polarons, resulting in a sudden collapse of the 
magnetic moment. Such a sharp drop signals the pres- 
ence of a localization mechanism beyond double exchange 
and is indeed seen experimentally (see, e.g., Refs. ||| puj ). 
On the other hand, at larger hole concentrations the po- 
laron binding energy is comparably small. Thus, a signif- 
icant suppression of the bandwidth via double exchange 
is needed before polaron formation can set in. The mag- 
netization curve now closely resembles the one predicted 
by double-exchange theory. 

Clearly beyond the grasp of conventional double- 
exchange theory lies the emergence of ferromagnetism 
in the insulating phase at low doping. Mostly responsi- 
ble for this are superexchange processes which mediate 
a ferromagnetic interaction even in the insulating phase. 
Ferromagnetism is further promoted by the existence of 
orbital polarons: Charge fluctuations inside the polaron 
provide a strong local ferromagnetic coupling between 
sites close to a hole, hence establishing ferromagnetic 
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clusters seen in experiment Jl3 At sufficiently large hole 
densities these clusters start to interact, thereby forming 
a ferromagnetic state. 

As was discussed above, orbital fluctuations are pre- 
dominantly induced by the motion of holes. The loss of 
charge mobility in the insulating phase should therefore 
trigger static orbital order. An orbitally ordered state 
has in fact been experimentally detected in the insulating 
regions of Lao.8sSro.12MnO3.E3 However, in general such 
an ordered state is expected to have orbital and Jahn- 
Teller glass features due to the presence of quenched or- 
bital polarons, thereby reducing the uniform component 
of Jahn- Teller distortions. Finally it is worth to notice 
that the phase diagram in this theory is highly sensitive 
to the transfer amplitude to as this parameter enters in 
the polaron binding energy. 



VI. CONCLUSION 

In summary, we have shown that a spontaneous devel- 
opment of orbital-lattice order is in general insufficient 
to trigger the localization process in manganites. Rather 
an additional mechanism was identified: Orbital polarons 
were illustrated to represent an intrinsic feature of an or- 
bitally degenerate Mott-Hubbard system and to play an 
important role in the physics of manganites. The binding 
energy of these orbital-hole bound states depends on the 
rate of orbital fluctuations and hence on the concentra- 
tion of doped holes: Polarons can form at low doping lev- 
els where orbitals fluctuate only weakly but they become 
unstable at higher levels of x. This scheme naturally in- 
troduces the hole concentration as an additional variable 
into the localization process. Most striking in this re- 
spect is the complete breakdown of metalicity observed 
below a critical hole concentration despite the fact that 
the system remains ferromagnetically ordered. On the 
other hand, orbital polarons become negligible at larger 
doping levels where the theory presented here converges 
onto a lattice-polaron double-exchange picture. Account- 
ing for both orbital and lattice effects we are finally able 
to reproduce well the important aspects of the phase di- 
agram of manganites. In general it can be concluded 
that a direct coupling between holes and surrounding or- 
bitals is of crucial importance for the physics of mangan- 
ites; its implications extend clearly beyond the metallic 
phase alone and can be expected to play an important 
role throughout the whole phase diagram. 
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